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ABSTRACT 

We have completed a 350 MHz drift scan survey using the Robert C. Byrd 
Green Bank Telescope with the goal of finding new radio pulsars, especially 
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millisecond pulsars that can be timed to high precision. This survey covered 
~10300 square degrees and all of the data have now been fully processed. We 
have discovered a total of 31 new pu l sars, seven of which are recycled pulsars. 



A companion paper by iBoyles et al.l ( 120121 ) describes the survey strategy, sky 



coverage, and instrumental set-up, and presents timing solutions for the first 
13 pulsars. Here we describe the data analysis pipeline, survey sensitivity, and 
follow-up observations of new pulsars, and present timing solutions for 10 other 
pulsars. We highlight several sources — two interesting nulling pulsars, an isolated 
millisecond pulsar with a measurement of proper motion, and a partially recycled 
pulsar, PSR J0348+0432, which has a white dwarf companion in a relativistic 
orbit. PSR J0348+0432 will enable unprecedented tests of theories of gravity. 

Subject headings: pulsars: individual (J0348+0432, J0458-0505, J1501-0046, 
J1518-0627, J1547-0944, J1853-0649, J1918-1052, J1923+2515, J2013-0649, 
J2033+0042)— surveys 



1. Introduction 

The vast majority of observed neutron stars in the Galaxy manifest themselves as radio 
pulsars. The extremely high rotational stability of pulsars, and especially millisecond pul- 
sars (MSPs), make them unrivaled laboratories for studying a wide range of astrophysical 
phenomena. Most pulsars have been discovered in large-area surveys, but most of these have 
focused on southern declinations or narrow regions around the Galactic plane. There is a 
need to find more pulsars in the northern sky, particularly high-precisi on MSPs that ca n 



be included in a pulsar timing array to detect gravitational waves (e.g. iJenet et al.ll2009l ). 
The 100-m Robert C. Byrd Green Bank Telescope (GBT) is one of the best telescopes in 
the world for finding and studying pulsars and a visible-sky pulsar survey using the GBT is 
underway (the Green Bank North Celestial Cap survey). 

During the northern summer of 2007 the azimuth track of the GBT underwent repair, 
making normal operations impossibld^l. Our team took advantage of the situation by com- 
pleting the GBT 350 MHz Drift-scan Survey between May and August. Because the GBT 
was unable to move in azimuth, we observed at a number of fixed elevations and allowed the 
sky to drift through the telescope beam at the sidereal rate. This survey was one of several 
low- frequency GBT surveys that are optimized for finding bright, nearby pulsars, with an 



15 A history of the track repair can be found at |http : //www . gb . nrao . edu/gbt/track . shtml 
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emphasis on MSPs. These surveys have either been completed (IHessels et al.l 120081 ) . are on- 
going, or are planned for the future. We collected over 1491 hours of data totaling 134 TB. 
Approximately 30 TB of these data are being analyzed by the Pulsar Search Collaborator}!^!, 
an educational initiative that actively involves high school stude nts and teachers in research 
under the guidance of a team of astronomers (IRosen et al.ll2010l ). For the remainder of this 
paper we discuss only the ~100 TB of the data that we have analyzed ourselves. All of 
these data have been fully processed and we have discovered 31 new pulsars, including ten 
recycled pulsars (seven of which are MSPs with P < 10 ms). We have derived full timing 
solutions for 2 5 of these new pul sars. The first 13 pulsars are being presented in a com- 
panion paper (IBoyles et al.ll2012l hereafter iPaper t ). along with a detailed description of 
the survey strategy, sky coverage, and instrumental set-up. We present timing solutions for 
an additional 10 pulsars here and describe the survey pipeline and data analysis in detail. 
An earlier-discovered Drift-s can pulsar, PSR J1023+0038, has been discussed elsewhere (see 
Archibald et~aD bood boioh . while another MSP, PSR J2256-1024, will be presented in a 
future paper (Stairs et al. in prep). In §|2]we explain how the data were divided into regions 
on the sky, interference removal, our de-dispersion scheme, and search strategies. In §|3]we 
describe our approximate survey sensitivity and the effects of scattering. In §4] we describe 
how we confirmed candidate pulsars and our follow-up observations. In §5] we present timing 
solutions and discuss some interesting individual systems. A summary can be found in £j6j 



2. Data Analysis 



The Drift-scan survey covered ~10300 square degrees. A deta iled des cription of the 
survey strategy, sky coverage, and instrumental set-up can be found in lPaper 1 Here we focus 
on the data re duction and se arch techniques. All data were processed using the PRESTC0 
software suite ( jRansomll200ll ). 



2.1. Pseudo-Pointings and Interference Excision 

Data were collected while the azimuth track of the GBT was being repaired, so the 
telescope was locked at constant azimuth. Different regions of the sky were observed by 
changing the elevation of the telescope and allowing the sky to drift through the telescope 



le 



http : //www . pulsarsearchcollaboratory . com/ 



1 ' http : //www . cv . nrao . edu/ -sransom/presto/ 
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beam at the sidereal rate. The time for a given point on the sky to pass through the beam 
is 

t^ - b (1) 
/<- s id cos 

where b = 36' is the full-width-at-half-maximum of the GBT beam at 350 MHz, 7£ S id = 
15' min -1 is the sidereal rate, and 5 is the declination. The survey covered —8° < 5 < 



+38° and —21° < 5 < +38°, depending on the azimuth of the telescope (see iPaper t for 
details). Although the telescope was not actually tracking the sky, we defined an individual 
pseudo-pointing to be a continuous block of data ~140 s in duration. Each pseudo-pointing 
overlapped with the preceding one by 70 s, so that all of our data were processed as part of 
two different pseudo-pointings. 



The raw data were collected using the GBT Pulsar Spigot back-end (IKaplan et al.ll2005l ). 
The Spigot uses autocorrelation chips that each work on 3-level raw samples and create 
an adjustable number of lags. These are then integrated into either 8- or 16-bit values, 
depending on the mode. The center frequency of the observations was 350 MHz and the 
bandwidth was 50 MHz. Most of our observations were made using the 8-bit mode that split 
the band into 2048 lags which were fast Fourier transformed to synthesize 2048 frequency 
channels, each with a width of 24.4 kHz, and recorded every 81.92 /is. A relatively small 
amount of data was taken early in the survey using the 16-bit, 1024 channel mode with the 
same sampling time. Hence, each 140 s block of data consisted of roughly 1.7 million spectra. 

Each pseudo-pointing was independently analyzed for radio frequency interference (RFI) 
using the rf if ind tool from PRESTO. Data were broken into blocks roughly two seconds long 
while maintaining the full frequency resolution. The total power, mean, and variance in both 
the time and frequency domain were calculated for each data block and compared to the 
median quantity for the entire pseudo-pointing. A time/frequency block was flagged as RFI 
and masked out (i.e. set to zero) in future analysis if the value of total power, mean power, 
or variance was greater than ten/four standard deviations from the median values for the 
entire pseudo-pointing. If greater than 30 %/70% of time intervals/channels were flagged then 
all remaining intervals/channels were masked out as well, under the assumption that they 
probably contained RFI just below our cut-off threshold. In addition to blindly searching for 
RFI, the Fourier spectra were de-reddened and persistent, well known sources of interference 
(e.g. the 60 Hz signal from AC power sources) were explicitly removed from the power 
spectra of each pseudo-pointing. Despite the fact that there was significant construction on- 
site due to the track repair, our data were remarkably free of RFI. For example, the median 
masking fraction was 0.56% and only 0.2% of the data had a masking fraction greater than 
30%. We are confident that we reached the noise limit for the vast majority of our survey 
(see §3). 
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2.2. De-Dispersion 



Free electrons in the interstellar medium give rise to a frequency- dependent dispersive 
time delay which, if left uncorrected, will make it virtually impossible to find new pulsars. 
The magnitude of the delay between two frequencies, i>\ and u 2 , is 



£ DM ^ 4.15 x 10 3 s x DM x 



MHz 



MHz/ 



(2) 



where DM is the dispersion measure in units of pc cm -3 . After applying appropriate shifts to 
each frequency channel, we summed over frequency to create de-dispersed time series. Each 
time series was transformed to the Solar System barycenter using the DE200 ephemeris 
(the default used by PRESTO) and the dispersion delay was removed (i.e., as if the signal 
had infinite frequency). We note that we used the DE405 ephemeris for baycentering when 
deriving pulse times of arrival (see §4.1(1 . 



The finite size of a frequency channel will still induce a smearing given by 



^chan 



.3 x 10 3 s 



/ DMAz/ chan \ 



(3) 



where Az/ C h an is the channel width and v is the channel center frequency (both in megahertz). 
For our primary observing mode with 24.4 kHz channels centered between 325 and 350 MHz, 
tchan ~ (3.8 to 5.9) fis x DM. 

Since the DM of a pulsar is not known a "priori, we created de-dispersed time series for 
each pseudo-pointing over a range of DMs, from to ~1000 pc cm -3 , which is a factor of 



three to four larger than the maximum DMs predicted by the NE2001 model (jCordes fc Lazio 
20021 ) in the low Galactic latitude regions of the survey. The step size between subsequent 
trial DMs (ADM) was chosen such that over the entire band ^adm — ^chan- This ensures 
that the maximum extra smearing caused by any trial DM deviating from the source DM 
by ADM is less than the intra-channel smearing 18 !. To increase computational efficiency, the 
data were down-sampled in time by adding 2 n samples together (where n is an integer) when 
Wn > 2 n x 81.92 /is. 



18 Dc-dispersion plans were generated using the DDplan.py tool in PRESTO. 
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2.3. Search Algorithms 

2.3.1. Periodic Sources 

Each RFI-cleaned, de-dispersed time series was Fourier transformed and searched for 
periodic signals. As mentioned in §2.1[ known sources of RFI were explicitly removed from 
the Fourier spectrum, so there is a very small chance (~0.006%) that a pulsar with a spin 
frequency very close to known RFI could have also been removed. Acceleration searches for 



isolated and binary pulsars were carried out in the Fourier domain ( IRansom et al.ll2002l ) for 
signals with a maximum drift of z max = ±50/nharm Fourier bins, where n harm is the highest 
harmonic where the pulsar is detected. This corresponds to a physical acceleration of 

z cP 

A - Zmax °- r (A) 
n max ,2 5 \*J 

where c is the spe ed of light, P is the sp in period of the pulsar, and t int = 140 s is the effective 
integration time (IRansom et al. H2OO2L For a P = 2 ms pulsar detected with up to eight 



harmonics, A max 24 ms~ 2 . Acceleration searches used up to eight summed harmonics, 
but we also carried out searches for unaccelerated pulsars (z = 0) using up to 16 summed 
harmonics. Only powers-of-two numbers of harmonics were summed. 

To filter out spurious and low-significance signals, only candidates that appeared in 
at least two time series of different DMs passed to the next stage of consideration. We 
also filtered duplicate signals (keeping only those with the highest S/N) that were within 
±1.1 Fourier bins of each other in different DM time series, as well as those that were 
harmonically related to each other. We folded up to 20 of the remaining candidates from the 
zero-acceleration searches and up to ten from the high-acceleration searches if their Fourier 
power was at least 6a above the Gaussian-equivalent noise level. We used the prepfold 
routine in PRESTO to fold the full resolution data at the nominal P, period derivative (P), 
and DM as determined by our searches. Our folding code refined these values and created 
diagnostic plots that were then saved for human inspection. 



2.3.2. Single-pulse Sources 

We searched for bright single pulses using single_pulse_search.py in PRESTO. Each 
time series was smoothed using a piecewise linear fit to the data, where each piece was 
2000 points long. The smoothed data were then correlated with boxcar functions of varying 
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widtha^j, which acted as matched filters to individual pulses. We recorded all single pulses 
with a signal-to-noise ratio, S/N > 5 and created diagnostic plots for all pulses with S/N > 
5.5. These plots were then saved for human inspection. In addition, an automated algorithm 
was used to flag pseudo-pointings with promising candidates, which were then inspected in 
greater detail. Five rotating radio transients (RRATs) have been discovered in this survey 
and a further 26 candidates have been identified and await confirmation. These discoveries 
and the automated algorithm used to help identify them will be presented in a forthcoming 
paper (Karako-Argaman et al. in prep). 



3. Survey Sensitivity 



Following iDewey et al.l (119851 ) and lLorimer fc Kramer! (120051 ) , the sensitivity of a pulsar 



survey may be written in terms of the phase-averaged limiting flux density 



o g (S/N mi jT sys / W 

where e is a degradation factor (discussed below), /3 = 1.16 is a correction factor that accounts 
for digitization losses, S/N min is the S/N threshold, T sys is the total system temperature, G 
is the telescope gain, n p is the number of summed polarizations, t- mt is the integration time, 
Av is the bandwidth, and W is the total pulse width (see Tabled] for relevant values). The 
degradation factor e accounts for drift of the pulsar through the telescope beam, which is 
not uniform in sensitivity. For a Gaussian primary beam 

eoc / e-^W^dt, (6) 
Jo 

where r(t) is the distance from the beam center and / = &/(2vTn 2). From simple geometry, 
r 2 (t) = y 2 + (b/2 — xt) 2 , where y and x are the distances from the beam center in right 
ascension and declination, respectively, and x is the drift rate. We normalize e such that a 
pulsar at the center of the beam for an entire integration has e = 1. For reference, a pulsar 
that crosses the beam center will have e = 0.81. 

The system temperature is a sum of several factors, including the receiver temperature 
(T rec ) and the sky temperature (T s k y ). The 350 MHz receiver^ of the GBT has a nominal 



19 The boxcar functions had a maximum width of either 150 x n x dt or 0.1 s, whichever was greater, where 
n is the down-sampling factor. 

20 We have folded the receiver spillover and cosmic microwave background into this number. Characteristics 
of the GBT receivers are available on the NRAO website (http://www.gb.nrao.edu/astronomers.shtml). 
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T r ec = 23 K. The Galactic synchrotron emission contributes heavily to T s k y , but this depends 
on sky position. Most of our survey w as at high Galactic latitudes, where the synchrotron 
emission adds ~30-50 K at 350 MHz (IHaslam et al.lll982l ). Our line of sight through the 
Galaxy also affects our sensitivity by increasing scattering and dispersion, both of which 
contribute to the observed pulse width. The typical maximum predicted DM at the high 
Galactic latitudes we cover is ~60 pc cm" 3 . According to the NE2001 model, this corre- 
sponds to a scattering time ~0.08 ms at 350 MHz, though observed scattering times may 
differ from predictions by an order of magnitude or more. Obviously, DM effects become 
much worse at low Galactic latitudes. Figure [1] shows approximate sensitivity curves for 
various combinations of y (the minimum offset from the beam center), T sys , and DM. These 
calculations do not take the effects of RFI into account, but as we describe in §2.1, the survey 
did not suffer greatly from RFI contamination. 



4. Candidate Confirmation and Follow-up 

Periodic and single-pulse candidates from each pseudo-pointing were judged by eye. 
Folded candidates were usually judged on three main criteria: i) distinct peaking of the 
signal's significance at DMs greater than pc cm -3 ; ii) broad-band emission (allowing for the 
possibility of regions of enhanced/diminished flux due to interstellar scintillation); and Hi) 
fairly persistent emission in time (allowing for eclipses and nulls and accounting for the roll- 
off in sensitivity near the edges of the telescope beam). In the case of single-pulse candidates, 
we looked for pulses that peaked at a non-zero DM and that decreased in significance away 
from this peak. Multiple single pulses at the same DM were also an obvious indicator of a 
good candidate. 

Promising candidates were confirmed in follow-up observations with the GBT, after 
which we began regular timing observations. To improve the quality of initial timing solu- 
tions, new pulsars had their sky positions refined by observing at a grid of locations with 



smaller GBT beams at successively higher frequencies ( IMorris et al.ll2002f ). We used a num- 
ber of dense observations early in the campaigns to characterize the orbits of binary pulsars. 
The majority of timing observations were carried out at 820 MHz, but most pulsars were 
also observed at other frequencies, allowing us to explore their spectral pro per ties. We 
also started using th e new Green Bank Ultimate Pulsar Processing Instrument! 21 ! (GUPPI) 



dDuPlain et all 120081 ) in 2008 October. Compared to the Spigot back-end, GUPPI offers 



larger bandwidth, better frequency and time resolution, higher dynamic range, and greater 



21- 



https : //safe . nrao . edu/wiki/bin/view/CICADA/NGNPP 
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resilience to strong RFI. 



4.1. Pulsar Timing Analysis 

All of the new pulsars were observed regularly with the GBT as part of our timing 
campaign. Long period pulsars (with P > 0.1 s) were observed for a minimum of about 11 
months, while MSPs and recycled pulsars were observed for a minimum of 20. Each pulsar 
was typically observed for 10-15 min per observing session. High-S/N average pulse profiles 
were created for each observing frequency by summing data from multiple observations. 
We created standard pulse profiles by fitting one or more Gaussians to these average pulse 
profiles using a least squares minimization routine^ 2 ! These standard profil es were used 



to co mpute pulse times of arrival (TOAs) using either PRESTO or PSRCHIVH 23 ! (IHotan et al. 



20041 ) (depending on data format) by cross-correlation in the Fourier domain. We typically 
obtained two TOAs per observation for isolated pulsars and four to six TOAs per observation 
for binary pulsars, ensuring goo d sampling of the orbit. Phase connected timing solutions 
were created using the TEMPO 24 ! software package and the DE405 Solar System ephemeris. 
All of our timing solutions are referenced to UTC(NIST). All the pulsars timed here have 
timing solutions with reduced x 2 > 1. Since we observe no unmodeled trends in our timing 
residuals, this is probably due to an underestimate of individual TOA uncertainties. As is 
standard practice, we multiplied all TOA uncertainties by an "error factor", such that the 
reduced x 2 — 1- 



4.2. Flux Measurements 

Mean flux densities were estimated by assuming that the off-pulse root-mean-square 
(RMS) noise level was described by the radiometer equation, 

a = f Tsys , (7) 

Gyra p Az/t int 

where T sys is the total system temperature. To ensure a proper estimate of the RMS noise 
level, we fit a third order polynomial to the off-pulse region and then subtracted this to 



pygaussf it .py in PRESTO 

http : //psrchive . sourcef orge . net/ 



24- 



http : //tempo . sourcef orge . net 
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create a flat off-pulse baseline. It is important to keep in mind, however, that observed 
pulsar fluxes are variable due to interstellar scintillation. The values that we report here are 
determined by averaging several observations but should be treated only as representative. 
We also calculated the spectral index when flux density estimates were available for multiple 
frequencies. We did this by fitting a standard power law to the flux density estimates, 
typically at 350 MHz and 820 MHz, assuming S u oc u a . The average for the pu lsars presented 



here i s (a) = —1.7, which is very similar to the average value presented in lLorimer et al. 



(119951 ) . 



We also attempted to measure the rotation measure (RM) whenever fully calibrated 
polarization data were available (which was at least once for each pulsar). We searched 
over a wide range of RMs, from ±1000 rad m -2 . We could only detect a significant RM 
for a subset of pulsars. Those pulsars without reported RMs are probably weakly polarized 
sources. 



Results 



A total of 31 n ew pulsa rs have been discovered thus far in the Drift-scan Survey. The first 
13 are presented in IPaper H and 10 are discussed here. A s ment ioned in jCD PSR J102 3+0038 
has been discussed elsewhere by (lArchibald et al.ll2009f ) and (lArchibald et al.ll2010l ). while 
another MSP, PSR J2256— 1024, will be presented in a future paper (Stairs et al. in prep). 
Full timing solutions have not been obtained for the six most recently discovered pulsars 
and these will be presented in future work. The 10 pulsars presented here include eight long 
period, isolated pulsars, one mildly recycled binary pulsar, and one isolated MSP. Seven of 
these 10 pulsars were detected independently in our searches for single pulses. Full timing 
solutions and other properties for the long-period pulsars are presented in Table |2] and for 
the recycled pulsars in Table [3j Integrated pulse profiles can be seen in Figures |2] and [3] 
and post-fit timing residuals are shown in Figure HJ We discuss some individual systems in 
greater detail below. 



5.1. PSR J0348+0432: A Relativistic Binary with a Low-Mass Companion 



PSR J0348+0432 (hereafter J0348) is a mildly recycled binary pulsar with P = 39.1 ms. 
The low magnetic field of J0348 (B gm { = 3.1 x 10 9 G) indicates that it is indeed partially 



- 11 - 



recycled and not a young pulsail 25 !. The DM of J0348 is 40.5 pc cm -3 and the DM-derived 
distance is 2.1 kpc. The orbital period of this system is 2.4 hr, and only three pulsars 
with P < 40 ms outside of globular clusters have shorter periods. If we assume a mass of 
1.4 M for J0348 then the observed mass function, f(M) = 2.9 x 10~ 4 , implies a minimum 
companion mass of 0.086 M . We searched for and identified an optical counterpart to 
J0348 in the Sloan Digital Sky Survey with corrected SDSS magnitudes u' = 21.84 ± 0.19, 
g' = 20.71 ± 0.03, r' = 20.60 ± 0.03, i' = 20.69 ± 0.05, and z' = 20.40 ± 0.15. More detailed 
spectroscopic follow-up with the Apache Point 3.5-m telescope and the Very Large Telescope 
have shown that the companion to J0348 is a low-mass white dwarf. The combination of a 
neutron star and low-mass white dwarf in a tight, relativistic orbit is unique among pulsars 
and makes J0348 an excellent laboratory for testing general relativity and other theories of 
gravity. Specifically, theories that invoke a scalar gravitational field predict that J0348 will be 
a strong emitter of dipolar gravitational radiation because the very different binding energies 
of th e neutron st ar and white dwarf will lead them to couple differently to the scalar fiel d 



Stairsl l2003l). Similar te sts have been done with PSR J1141-6545 fahat et al.l l2008h . 



(e.g. 

J1012+5307 (lLazaridis et al.l l2009h and J1738+0333 flFreire et al.ll2012l ). but J0348 is in a 
tighter, more relativistic orbit and likely has a less massive companion, so it will be a stronger 
probe of these effects. A full analysis of the spectroscopic observations of J0348 and their 
implications for alternative theories of gravity will be presented in Antoniadis et al. (in 
prep). 

J0348 shows significant profile evolution as a function of frequency (see Figure [3]); at 
frequencies above 1.4 GHz, the main profile component becomes extremely narrow, with 
a duty cycle of ~1%. This has allowed us to obtain very precise pulse arrival times — our 
RMS timing residuals are 9.3 /is but at high frequencies individual TOA uncertainties can 
be < 3 /is. We are continuing long-term timing of this pulsar using the Arecibo Observatory. 



5.2. PSR J1923+2515 

PSR J1923+2515 (hereafter J1923) is an isolated MSP with a 3.8-ms spin period. Fig- 
ure [3] shows the integrated pulse profile of J1923 at several different frequencies and the 
evolution in the profile shape is clear. We see evidence for a weak interpulse in our summed 
820 MHz and 2 GHz data. The timing of J1923 improved significantly at higher frequencies, 
where we were able to obtain TOAs with uncertainties < 1 /is. J 1923 is being regularly 



25 Although we have not measured proper motion and cannot calculate its contribution to P hs, it is 
certainly not sufficient to increase the i? sur f by several orders of magnitude. 
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obser yed at Arecibo as pa rt of the NANOGrav timing array for gravitational wave detec- 



tion ( D cmorestetaL 



( Ivan Haasteren et al. 



2012h . It will also be suitable for the European Pulsar Timing Array 



2011) 



J1923 is the only pulsar presented here for which we were able to measure a significant 



proper motion. We find /i c 



-6.2(1.2) masyr 1 and n$ = —23.6(3.5) masyr 



We 

performed an F-test to determine if the addition of proper motion is in fact required by the 
data. The full x 2 °f our timing model excluding proper motion is 339.63, with 147 degrees of 
freedom. When proper motion is included in the fit, x 2 — 146.02 with 145 degrees of freedom. 
The probability that this improvement is due to chance is 1 x 10~ 16 . Thus, the improvement 
in our timing solution when proper motion is included is extremely significant. We used 
the DM and the NE2001 model of Galactic free electron density to estimate the distance to 
J1923, D = 1.6(3) kpc, where the nu mber in parentheses r epresents a 20% fractional error, 
which is typical for these estimates (ICordes fc Laziol 120021 ) . At this distance, the observed 
proper motion corresponds to a transverse velocity v± = 188(4 6) kms" 1 , which is within 



the observed range of other MSPs, though higher than average (IGothoskar fc Gupta! 12000 
Boedanov et al.ll2002[ bonzalez et aDboilh . 



Using this v± we can calculate the magnitude of the Shklovskii effect (1 S hklovskiil 1 1 9 701 ) : 



P 



JL 
P 



cD 



We find P M = 8.9(4.7) x 10 21 s s 1 . Acceleration within the Galac tic potential wi l l also cause 



a bias in the observed P. We estimate this contribution following iNice fc Taylor! ( 119951 ). but 
find that biases due to acceleration perpendicular and parallel to the Galactic plane are only 
—3.3 x 10 -22 ss _1 and —3.4 x 10~ 22 ss" 1 , respectively. These are an order of magnitude 
smaller than P^. The bias due to the Shklovskii effect is 94% of the observed P and would 
imply that the intrinsic spin-down of the pulsar is P; nt = 0.6(4.7) x 10 -21 ss -1 , so we can 
only place an upper limit on P int at this time. A more precise measurement of proper motion 
or a better distance estimate will be needed to constrain the magnitude of the Shklovskii 
effect and to obtain a better measurement of P m t. In the meantime, the derived quantities 
listed in Table [3] are upper limits based upon our measurement uncertainties for Pj nt . 



5.3. PSR J0458-0505 

PSR J0458— 0505 (hereafter J0458) is a nulling pulsar with a 1.9-s spin period. It 
was detected in both the Fourier domain and single-pulse searches. The profile is slightly 
asymmetric (see Figure [2]), with a small trailing component. The reduced x 2 obtained from 
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our timing solution after fitting for position, P, P, and DM was substantially higher than 
unity. As we see no systematic trends in our timing residuals, we assume that the number 
of pulses per individual observation was too small, due to a combination of a long period, 
limited integration time and a very large nulling fraction (see below). The pulse profile thus 
probably did not stabilize within these TOA integrations. To make the reduced \ 2 equal 
to one, we multiplied our individual TOA errors by a constant factor of 2.97. Although we 
were still able to derive an accurate timing solution for J0458, the fractional errors in the 
timing parameters are larger than for most of the other pulsars presented here, especially 
for declination and P. 



5.3.1. Estimate of the Nulling Fraction 

We estimate the nulling fraction, NF, of J0458 in the following way. We first removed 
strong sources of RFI from each observation. We then folded each data-set using sub- 
integrations that were a single pulse period in duration using the psrf its_singlepulse 
command from psrf its_utila 26 l. In some sub- integrations, systematic trends due to lower 
levels of RFI were still visible. To remove these, we used a least-squares minimization 
technique to fit up to a maximum of four independent sinusoids to the off-pulse region 
of each sub-integration and then subtracted them from the data, creating a flat off-pulse 
baseline. Each sub-integration was then normalized to have an off-pulse median and RMS 
noise level of zero and one, respectively. 

To determine if the pulsar was in an "on" state, we calculated the integrated S/N in the 
on-pulse region, which was determined by inspection of the integrated pulse profile. We also 
calculated the integrated S/N in an off-pulse region with the same number of bins as the 
on-pulse region. We counted a pulse as being in the "on" state if it had an integrated S/N 
above some threshold, S/N thresh . We chose S/N thresh based on the statistics of the off-pulse 
region. Histograms of the on-pulse and off-pulse S/Ns can be found in Figure [5] and in Figure 
Owe show NF as a function of threshold S/N. Our calculations show that the off-pulse region 
rarely exceeded S/N = 3.0, as expected for Gaussian distributed noise. To be conservative, 
we set S/N thresh = 3.5, though we also report NF for S/N thresh = 3.0 and 4.0 for comparison. 

After removing RFI, J0458 was observed for a total of 2218 full rotations at 820 MHz 
and we find NF = [0.60,0.63,0.66] for S/N thresh = [3.0,3.5,4.0], respectively. J0458 was 
observed for a total of 978 rotations at 350 MHz with NF = [0.66, 0.69, 0.73] for S/N thresh = 
[3.0,3.5,4.0], respectively. These NFs are fairly high compared to other pulsars, but are not 



2e https : //github . com/demorest/psrf its_utils 



-14- 



unprecedented — PSRs B1112+50 and B1944+17 have comparable NF (jRitchingd Il976l ). 
The nulling fraction of J0458 seems to be similar at both 350 MHz and 820 MHz. This 
behavior is consistent with previou s studies that suggest nulling is a broadband phenomenon 
at low frequencies (see iBiggsl Il992l . and references therein) . 



5.4. PSR J2033+0042 



PSR J2033+0 042 (hereafter J2033) is a 5.0-s pulsar that, like J0458, nulls significantly. 
It was reported by iBurke-Spolaor Bailed (12 0101) as a RRAT that sometimes was detected 
in Fourier searches. IBurke-Spolaor fc Bailed ( 120101 ) report on the position, period, and DM 
of J2033 and note its high nulling fraction and the presence of drifting sub-pulses. Here, 
we present a full timing solution and quantitative measurement of the nulling fraction. We 
detected J2033 in both single-pulse and Fourier searches. Only nine radio pulsars and 12 
RRATs in the ATNF catalog have a longer period than J2033. Like J0458, the fractional 
errors in our timing parameters were relatively large, probably because our integration times 
were shorter than the pulse stabilization time. J2033 has a longer period than J0458 and 
also nulls significantly. 



We used the same procedure as outlined in §5.3.11 to estimate the nulling fraction 
in J2033. We observed the pulsar for 994 rotations at 820 MHz and we find a NF = 
[0.53,0.56,0.58] for S/N thrcsh = [3.0,3.5,4.0], respectively. J2033 was observed for 111 rota- 
tions at 350 MHz and we find NF = [0.44,0.48,0.49] for S/N thrcsh = [3.0,3.5,4.0], respec- 
tively. Histograms of S/N can be found in Figure [7] and NF as a function of S/N thresh is 
plotted in Figure [SJ The nulling fraction of J2033 is somewhat higher at 820 MHz than 
at 350 MHz, but we only observed J2033 at 350 MHz on one occasion, so a more detailed 
study of the nulling characteristics of this pulsar should be conducted before drawing firm 
conclusions about the frequency dependence of NF. Overall, pulsars J0458 and J2033 are 
on l ess than half of th e time, adding to a growing sub-population of pulsars that are mostly 
off faeane et alJboilh . 



6. Conclusion 

The Drift-scan Survey has discovered 31 pulsars, of which 10 are presented here. The 
majority are isolated long-period pulsars. J0348 is a mildly recycled binary pulsar that has 
a low-mass white dwarf companion in a relativistic orbit. It is a unique and powerful system 
for testing gravitational theories and hence we are continuing to time it long-term. A more 
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detailed study of J0348 will be presented in Antoniadis et al. (in prep). J1923 is an isolated 
MSP. We have a significant measurement of the pulsar's proper motion, but the implied 
magnitude of the Shklovskii effect is nearly equal to the observed spin-down, so we are only 
able to set limits on the rotational characteristics of J1923. Long-term monitoring should 
help to better constrain the proper motion and intrinsic spin-down. J0458 and J2033 are 
both nulling pulsars with nulling fractions > 50%. 
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Fig. 1. — Approximate phase-averaged limiting flux density of our survey. Black curves 
are for T sys = 75 K and DM = 30 pc cm~ 3 , while gray curves are for T sys = 100 K and 
DM = 75 pc cm -3 . The smallest offset from the center of the telescope beam is y. 
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Table 1. Parameters of the GBT 350 MHz Drift-scan Survey 



Paramter 


Value 


ADC Conversion Factor, /3 


1.16 


Signal-to- noise threshhold, S/N min 


6.0 


Receiver Temperature, T rcc (K) 


46 


Telescope Gain, G (K Jy" 1 ) 


2.0 


Number of summed polarizations, n p 


2 


Length of pseudo-pointing, t- m t (s) 


140 


Bandwidth, A/ (MHz) 


50 


Number of frequency channels, n c han 


2048 a 


Sampling time, i samp (ps) 


81.92 



a A small amount of data was recorded with 
1024 frequency channels early in the survey. 
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Table 2. Parameters of Newly Discovered Long Period Pulsars 



Parameter PSR J0458-0505 PSR J1501-0046 PSR J1518-0627 



Timing Parameters 



Right Ascension (J2000) 


04:58:37.121(26) 


15:01:44.9558(94) 


15:18:59.1104(80) 


Declination (J2000) 


-05:05:05.1(4.0) 


-00:46:23.52(88) 


-06:27:07.70(66) 


Pulsar Period (s) 


1.88347965849(18 ) 


0.4640368139284(82) 


0.7949966745699(78) 


Period Derivative (ss _1 ) 


5.3(1.5)xl0- 16 


2.391(60)xl0~ 16 


4.179(56)xl0~ 16 


Dispersion Measure (pc cm~ 3 ) 


47.806(32) 


22.2584(90) 


27.9631(98) 


Reference Epoch (MJD) 


55178.0 


55170.0 


55170.0 


Span of Timing Data (MJD) 


55006-55349 


55006-55335 


55006-55335 


Number of TOAs 


22 


39 


43 


RMS Residual (/is) 


716 


172 


175 


Error Factor 


2.967 


1.095 


1.050 


Derived Parameters 


Galactic Longitude (deg) 


204.14 


356.58 


355.15 


Galactic Latitude (deg) 


-27.35 


48.05 


41.0 


DM Derived Distance (kpc) 


2.6 


1.4 


1.6 


Surface Magnetic Field (10 12 Gauss) 
Spin-down Luminosity (10 32 ergs -1 ) 


1.0 


0.33 


0.58 


0.032 


0.95 


0.33 


Characteristic Age (Myr) 


56 


31 


30 


820 MHz FWHM 


0.014 


0.022 


0.012 


820 MHz Flux Density (mJy) 


0.5 


0.3 


0.4 


Spectral Index 


-1.6 


-2.1 


-1.8 
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Table 2. Parameters of Newly Discovered Long Period Pulsars (Cont.) 



Parameter PSR J1547-0944 PSR J1853-0649 PSR J1918-1052 



Timing Parameters 



Right Ascension (J2000) 


15:47:46.058(36) 


18:53:25.422(36) 


19:18:48.247(13) 


Declination (J2000) 


-09:44:7.8(3.2) 


-06:49:25.9(2.6) 


-10:52:46.38(66) 


Pulsar Period (s) 


1.576924632943(44) 


1.048132105087(54) 


0.798692542358(15) 


Period Derivative (s s _1 ) 


2.938(36) xlO" 15 


1.548(44) xlO" 15 


8.653(15) xl0~ 16 


Dispersion Measure (pc cm~ 3 ) 


37.416(22) 


44.541(36) 


62.73(80) 


Reference Epoch (MJD) 


55170.0 


55170.0 


55026.0 


Span of Timing Data (MJD) 


55006-55335 


54976-55335 


54712-55339 


Number of TOAs 


22 


24 


26 


RMS Residual (fjs) 


338 


588 


394 


Error Factor 


1.630 


2.150 


2.480 


Derived Parameters 


Galactic Longitude (deg) 


358.31 


27.08 


26.23 


Galactic Latitude (deg) 


33.57 


-3.55 


-10.96 


DM Derived Distance (kpc) 


1.9 


1.5 


2.1 


Surface Magnetic Field (10 12 Gauss) 


2.2 


1.3 


0.84 


Spin-down Luminosity (10 32 ergs" 1 ) 


0.30 


0.53 


0.67 


Characteristic Age (Myr) 


8.5 


11 


15 


820 MHz FWHM 


0.019 


0.015 


0.015 


820 MHz Flux Density (mJy) 


0.4 


0.5 


0.4 


Spectral Index 


-1.8 


-2.3 




Rotation Measure (radm~ 2 ) 




34.7(3.2) 


-47.6(6.0) 
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Table 2. Parameters of Newly Discovered Long Period Pulsars (Cont. 



Parameter 



PSR J2013-0649 PSR J2033+0042 



Timing Parameters 



Right Ascension (J2000) 
Declination (J2000) .... 



20:13:17.7507(38) 
-06:49:05.39(32) 



Pulsar Period (s) 0.5801872690010(34) 



Period Derivative (ss^ 1 ) 

Dispersion Measure (pc cm~ 3 ) 

Reference Epoch (MJD) 

Span of Timing Data (MJD) . 

Number of TOAs 

RMS Residual (jjta) 

Error Factor 



6.007(24) x 10- 
63.36(10) 
55172.0 
55005-55339 
45 
150 
1.370 



-14 



20:33:31.11(12) 
00:42:22.0(8.0) 
5.01339800063(90) 
1.013(78) x 10- 
37.84(13) 
55172.0 
55005-55339 
22 
2195 
4.680 



Derived Parameters 



Galactic Longitude (deg) 


36.17 


45.88 


Galactic Latitude (deg) 


-21.29 


-22.2 


DM Derived Distance (kpc) 


3.0 


1.9 


Surface Magnetic Field (10 12 Gauss) 


0.60 


7.2 


Spin-down Luminosity (10 32 ergs -1 ) 


1.2 


0.032 


Characteristic Age (Myr) 


15 


7.8 


820 MHz FWHM 


0.017 


0.018 


820 MHz Flux Density (mJy) 


0.6 


1.2 


Spectral Index 




-1.7 


Rotation Measure (radm~ 2 ) 




-71.2(2.2) 



Note. — Numbers in parentheses are 1-tr uncertainties as determined by TEMPO; 
although we have scaled the TOA uncertainties by the error factors reported, we have 
not doubled the nominal TEMPO uncertainties as is sometimes done in these cases. Flux 
density estimates typically have a 20-30% relative uncertainty due to scintillation. All 
timing solutions use the DE405 Solar System Ephemeris and the UTC(NIST) time 
syst em. Derived quantities as sume an R — 10 km neutron star with / = 10 45 gmcm 2 
(see lLorimer h Kramerl 120051 1 . The DM derived distances were calculated using the 
NE2001 model of Gal actic free electron density, and have typical errors of ~ 20% 
( Cordes fc Lazioll2002h . 
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Tabic 3. Parameters of Newly Discovered Short Period Pulsars 



Parameter 


PSR J0348+0432 


PSR J1923+2515 




Timing Parameters 


Right Ascension (RA; J2000) . . . 


03:48:43.63817(33) 


19:23:22.494560(76) 


Declination (DEC; J2000) 


04:32:11.449(10) 


25:15:40.6436(14) 


RA Proper Motion (masyr -1 ) . . 




-6.2(2.4) 


DEC Proper Motion (masyr -1 ) . 




-23.5(7.0) 


Pulsar Period (s) 


0.039122656280156(10) 


0.00378815551961303(52) 


Period Derivative (s s _1 ) 


2.417(16)xl0~ 19 


9.42(14) xl0~ 21 


Dispersion Measure (pc cm~ 3 ) . . 


40.56(11) 


18.85766(19) 


Reference Epoch (MJD) 


55278.0 


55322.0 


Span of Timing Data (MJD) 


54873-55682 


55005-55639 


Number of TOAs 


183 


153 


Error Factor 


1.657 


1.245 


RMS Residual (/js) 


10.33 


5.0 


Binary Parameters 


Binary Model 


ELL1 




Orbital Period (days) 


0.10242406134(30) 




Projected Semi-major Axis (lt-s) 


0.1409842(34) 




Epoch of Ascending Node (MJD) 


54889.70532337(65) 




1 st Laplace Parameter 


< 5.0 x 10~ 5 




2 nd Laplace Parameter 


< 6.3 x 10~ 5 




Derived Parameters 


Orbital Eccentricity 


< 8.1 x 10~ 5 




Mass Function (M ) 


0.000286807(20) 




Minimum Companion Mass (M Q ) 


0.086 




Galactic Longitude (deg) 


183.34 


58.95 


Galactic Latitude (deg) 


-36.77 


4.75 


DM Derived Distance (kpc) 


2.1 


1.6 


Transverse Velocity (km s _1 ) 




188(46) 


Shklovskii Effect (ss^ 1 ) 




8.9(4.7)xl0- 21 


Intrinsic Spin-down (ss _1 ) 




< 5.3 x 10~ 21a 
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Table 3 — Continued 



Parameter 


PSR J0348+0432 


PSR J1923+2515 


Surface Magnetic Field (10 9 Gauss) 


3.1 


< 1.4 a 


Spin-down Luminosity (10 32 ergs^ 1 ) 


1.6 


< 38 a 


Characteristic Age (Gyr) 


2.6 


> ll a 


820 MHz FWHM 


0.016 


0.142 


820 MHz Flux Density (mJy) 


1.8 


0.6 


Spectral Index 


-1.2 


-1.7 


Rotation Measure (radm~ 2 ) 


49.5(13) 


10.8(3.8) 



a These quantities are limits based on our error for P; nt after correcting for the 
Shklovskii effect. See £15.21 for further details. 



Note. — Numbers in parentheses are l-er uncertainties as determined by TEMPO; 
although we have scaled the TOA uncertainties by the error factors reported, we 
have not doubled the nominal TEMPO uncertainties as is sometimes done in these 
cases. Flux density estimates typically have a 20-30% relative uncertainty due 
to scintillation. All timing solutions use the DE405 Solar System Ephemeris and 
the UTC(NIST) time system. Derived quantities assume an R = 10 km neutron 
star with / = 10 45 gmcm 2 (see lLorimer fc Kramerl 120051 ) . Minimum companion 
masses were calculated assuming a 1.4 M Q pulsar. The DM derived distances 
were calculated using the N E2001 model of Galac tic free electron density, and 
have typical errors of ~ 20% ([Cordes fc Lazioll2002l) . 
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J0458-0505 

1.883 s 
0.6 mJy 



J1501-0046 
0.464 s 
0.3 mJy 



J1518-0627 
0.795 s 
0.5 mJy 



J1918-1052 

0.799 s 
0.5 mJy 



J 1547-0944 

1.577 s 
0.4 mJy 



J2013-0649 

0.580 s 
0.8 mJy 



J1853-0649 
1.048 s 
0.6 mJy 



J2033+0042 

5.013 s 
1.4 mJy 



Fig. 2. — Integrated 820 MHz pulse profiles for the newly discovered long period pulsars 
presented here. All profiles show one full rotation of the pulsar (i.e., from phase 0-1) with 
256 phase bins. The profiles were made by adding all the RFI-free observations for each 
pulsar and were used to create standard pulse profiles at 820 MHz. Pulse periods and mean 
flux densities are also given. 
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Fig. 3. — Integrated pulse profiles at four observing frequencies for the newly discovered 
recycled pulsars presented here. All profiles show one full rotation of the pulsar (i.e., from 
phase 0-1) with 256 phase bins. The summed profiles were made by aligning each folded 
profile using the TEMPO ephemeris and then adding all the RFI-free observations at the 
specified frequencies. The profile evolution in both pulsars is clear. The bars show the 
relative timescale for dispersive smearing at each frequency (see Equation [3]). 
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Fig. 4. — Post-fit timing residuals for each of the newly discovered pulsars. Only phase 
connected TOAs are shown. Nominal TOA errors have been multiplied by a constant error 
factor so that the reduced \ 2 — 1- Note that the axes have different scales in most plots. 
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Fig. 5. — Histograms of integrated S/N for single pulses of J0458. The top panel is for data 
taken at 350 MHz, while the bottom panel is for data taken at 820 MHz. The on-pulse 
region is shown in gray and the off-pulse region as hatched. The dashed line shows the 
S/N thrcsh = 3.5, above which we counted the pulsar as being in an "on" state. 
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Fig. 6. — The nulling fraction of J0458 as a function of S/N thresh . The top panel is for data 
taken at 350 MHz, while the bottom panel is for data taken at 820 MHz. 
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Fig. 7. — Histograms of S/N for single pulses of J2033. The data labels are the same as in 
Figure El 
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Fig. 8. — The nulling fraction of J2033 as a function of S/N thresh . The data labels are the 
same as in Figure EJ 



